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Light beams having a vectorial field structure 
— or polarization — that varies over the trans¬ 
verse profile and a central optical singularity are 
called vector-vortex (VV) beams and may ex¬ 
hibit specific properties, such as focusing into 
“light needles” or rotation invariance, with appli¬ 
cations ranging from microscopy and light trap¬ 
ping to communication and metrology. Individual 
photons in such beams exhibit a form of single¬ 
particle quantum entanglement between different 
degrees of freedom. On the other hand, the quan¬ 
tum states of two photons can be also entangled 
with each other. Here, we combine these two con¬ 
cepts and demonstrate the generation of quan¬ 
tum entanglement between two photons that are 
both in VV states — a new form of quantum “en¬ 
tangled entanglement”. This result may lead to 
quantum-enhanced applications of VV beams as 
well as to quantum-information protocols fully ex¬ 
ploiting the vectorial features of light. 

Quantum entanglement [ 1 ] lies at the basis of funda¬ 
mental questions on the nature of reality, as exempli¬ 
fied by the Einstein-Podolsky-Rosen argument [2] or the 
Schrodinger’s cat paradox [3]. On the other hand entan¬ 
glement is today also a key tool of quantum information 
technology, enabling applications such as quantum tele¬ 
portation [4] , superdense coding [5] and quantum compu¬ 
tation [ 6 ]. Entangled quantum states are nowadays rou¬ 
tinely produced with different physical systems spanning 
from atoms [7] to crystals [ 8 ] and photons [9]. Entangled 
photon pairs, in particular, are commonly generated by 
exploiting nonlinear optical processes [ 10 ] and may show 
entanglement in several degrees of freedom such as fre¬ 
quency, path, orbital angular momentum (0AM), and 
polarization [ 11 ]. 

The optical polarization - defining the oscillation di¬ 
rections of the electromagnetic fields - is typically ap¬ 
proximately uniform in a light beam. Yet, the polar¬ 
ization can also vary over the transverse profile, giving 
rise to vector beams with peculiar polarization patterns 
[12]. The so-called vector vortex (VV) beams are a par¬ 
ticular class of vector beams characterized by a cen¬ 
tral optical singularity surrounded by an azimuthally- 
varying pattern of polarization [13, 14]. These beams can 
be conveniently described as balanced nonseparable su¬ 
perpositions of polarization-0AM eigenmodes [15], with 


the 0AM magnitude defining the “order” of the beams. 
Hence, photons in VV beams are actually entangled in 
these two degrees of freedom: such kind of single-particle 
entanglement is also known as “intrasystem” entangle¬ 
ment and exists in classical systems too [16, 17]. VV 
beams have already found applications in areas ranging 
from microscopy [18] to metrology [19, 20], optical trap¬ 
ping [ 21 ], nano-optics [ 22 ], and quantum communication 
[23-25]. Due to their interesting properties, in the last 
years several techniques have been developed to generate, 
manipulate and analyze VV beams [13-15, 26, 27]. 

In this paper, we report the generation and charac¬ 
terization of entangled pairs of VV photons of arbitrary 
order. In particular, we consider five combinations of 
VV mode orders, corresponding to different polarization 
patterns for the two beams. We simultaneously demon¬ 
strate, by a complete sixteen-dimensional quantum to¬ 
mography, both the intrasystem entanglement between 
polarization and 0AM within each photon and the inter¬ 
system entanglement between the two photon states: the 
former is related to the structure of VV states, the latter 
corresponds to entanglement between two complex vecto¬ 
rial fields that has, to our knowledge, never been reported 
before. In this way, we are demonstrating experimen¬ 
tally a novel form of “entangled entanglement” [28, 29], 
namely hybrid entangled entanglement corresponding to 
two (intersystem) entangled pairs of (intrasystem) entan¬ 
gled pairs of qubits. Einally, by performing a non-locality 
test directly in the VV space, we show that entangle¬ 
ment between complex vectorial fields can be effectively 
exploited as a resource in fundamental quantum men- 
chanics as well as quantum information. 

Let us denote with \R,ti (|L,Q) the state of a photon 
with uniform right (left) circular polarization carrying Ih 
of orbital angular momentum. A VV beam of order m 
is defined in the two dimensional Hilbert space spanned 
by {|i?, m), |I/, —m)}. In particular, we consider the two 
balanced superpositions |r^) = ^ {\R,m) + \L,—m)), 

\'dm) = {\R,ni) — \L, —m)). Eor m 7 ^ 0, all these 

modes feature a polarization singularity and a null inten¬ 
sity in the center and hence exhibit the so-called dough¬ 
nut beam profile. When analyzed through a linear po¬ 
larizer, the corresponding intensity pattern consists of 
2m petals [19, 30]. When m = 1, the states |ri) and 
I'd!), corresponding to the well-known radially and az- 


2 


imuthally polarized beams [13], are also invariant un¬ 
der azimuthal rotations [23, 24]. For brevity, we will 
refer to states |f^) and |'d^) with the terms “radial” 
and “azimuthal”, irrespective of m. A generic VV beam 
(or photon VV state) can be represented on a “hybrid 
Poincare sphere” (HPS) [31], also known as “higher-order 
Poincare sphere” [32], where states {\R,m ), jL, —m)} lie 
on the poles and |f^) and \'dm) he on opposite points 
on the equator. Figure 1 shows examples of this rep¬ 
resentation and the polarization and intensity patterns 
of some of these modes. By noticing that the complete 
polarization-0AM Hilbert space of order m is spanned 
by the four states {|i?, m), jL, —m), jL, m), |i?, —m)}, 
we can also define 7r-modes of order m as bal¬ 
anced superpositions in the Hilbert space spanned by 
{\R,-m) ,\L,m)}: |7r+) = {\R,-m) + \L,m)) and 

l^m) = 71 - \L,m)) 

In this work, we generate and detect VV modes by us¬ 
ing a specialized optical component, the g-plate, a bire- 
fringent patterned slab that can couple/decouple polar¬ 
ization and 0AM of single photons [33]. More in de¬ 
tail a g-plate with topological charge q maps a photon 
with input state a\R,0) /3 \L,0) into the output state 

a |I/, —2q) -h /3 |i?, 2q) and vice versa. Thus, radial and 
azimuthal VV beams |f^) and \'3rn) are easily produced 
by using a linear horizontal (H) and vertical (V) input 
polarization, respectively, with m = 2q. Superpositions 
of VV beams are also obtained from the corresponding 
input polarization superpositions. A g-plate can also be 
used to measure VV beams. Indeed, a radial (azimuthal) 
VV beam of order m is transformed into a linear hori¬ 
zontal (vertical) uniformly polarized beam by a g^-plate 
with 2q = m. In this way, the measurement of a com¬ 
plex polarization pattern such as that of a VV beam is 
reduced to a much simpler polarization measurement, be¬ 
cause the g'-plate acts as an “interface” between the two 
dimensional space of VV beams of order m and the po¬ 
larization one [14, 27]. A more general and complete ap¬ 
proach to analyze VV beams is by measuring separately 
the polarization and 0AM, for instance by performing a 
quantum state tomography in the complete polarization- 
0AM Hilbert space. 

The experimental apparatus we used for the generation 
of entangled photons in VV modes is depicted in Fig. 2 
(Generation section). A polarization-entangled photon 
pair is produced by exploiting spontaneous parametric 
down conversion in a /3-barium borate crystal (BBO) [10]. 
The generated state is (|^)i 1^)2 “ l^)i 1 ^) 2 )’ where 
the subscript denotes the photon output path. After the 
crystal, both photons are filtered in wavelength, via a 
narrowband interference filter, and in spatial mode, via a 
single-mode fiber. This last operation sets the transverse 
mode of both photons to a fundamental Gaussian one 
(TF/Mqo), corresponding to vanishing 0AM. As a last 
step, two q-plates with topological charges qi and g 2 , 



Figure I. Representation of entangled VV beams and their 
underlying hybrid entangled entanglement structure. Panel 
(A) shows the intersystem entanglement (green dashed lines) 
between VV beams belonging to different orders. Each VV 
state is defined on a hybrid-Poincare sphere: purple spheres 
HP Si and HPS 2 , corresponding to orders m — 1 and m = 5 
respectively, represented together with the intensity and po¬ 
larization patterns for the corresponding azimuthal, radial 
and circularly-polarized modes (for graphical reasons, these 
spheres are oriented with the two poles corresponding to 
polarization-0AM eigenstates aligned in the horizontal direc¬ 
tion). The Poincare spheres in the two lower boxes represent 
the polarization tt (blue spheres) and 0AM o\^ (red spheres) 
spaces within each photon, with their respective polarization 
patterns and phase profiles. The intrasystem entanglement 
(blue/red dashed lines) between these two degrees of freedom 
of the single photons generates the VV states. Panel (B) re¬ 
ports the intensity and polarization patterns of radial and 
azimuthal VV beams of various orders. 


transform the two polarization photon states into VV 
modes of orders mi = 2^^! and m 2 = 2 ^ 2 , respectively. 

The resulting state emerging from the source is the 
following: 

1^7711,7712) ~ |'^ 7772)2 ~ 1^7712)2) ’ (^) 

which corresponds to a pair of entangled vector-vortex 
beams of orders mi and m 2 . A complete set of VV Bell 
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Figure 2. Experimental apparatus for the generation and 
analysis of entangled VV photons. In the Generation stage, 
a polarization-entangled photon pair is produced in a BBO 
crystal and is then converted into a VV-entangled pair by two 
^-plates. The generation of a 7r-mode, needs an additional 
HWP after the g-plate. The Tomo section corresponds to 
a polarization-analysis stage formed by a quarter-wave plate 
(QWP), a HWP and a polarizer (PBS). This section is used 
only for the full tomography in the complete polarization- 
0AM space, while it is not needed for direct VV-space mea¬ 
surements. The Analysis stage deals with the 0AM- or VV- 
mode analysis (it is OAM-mode if the Tomo-stage is inserted, 
otherwise it is VV-mode): on each arm a ^-plate converts 
back the 0AM- or VV-modes into uniform polarization states, 
which are then measured with a polarization analysis stage. 
The photons are then coupled into a single mode fiber that 
filters the non-Gaussian modes and are sent to single photon 
detectors. 


states can be obtained by performing local operations on 
one of the photons. Indeed an input state in polar¬ 
ization is transformed by the two g-plates into 
and |^+) is transformed into same applies 

for any linear combination of Bell states: this versatile 
technique allows to generate arbitrary maximally entan¬ 
gled vector vortex beams. 

To prove the generality of our approach, we generated 
various combinations of VV mode orders. By exploiting 
g'-plates with q = 1/2, 3/2, 5/2 and 5 we considered the 
following combinations (mi,m 2 ): (1,1), (1,5) (1,10) and 
(3,5). Finally, we generated the pair (1,-1). The VV 
modes with m = — 1 are 7r-modes [34], which can be 
obtained by flipping the circular polarization handedness 
with a half-wave plate (HWP) added after the g-plate 
with q = 1/2, or by directly exploiting a g-plate with 
q = -l/2. 

To fully analyze the entangled pairs, we performed a 
quantum state tomography in the polarization and the 
0AM subspaces corresponding to VV mode of each pho¬ 
ton. In this way it is possible to measure both the 
intrasystem and the intersystem entanglement of our 
states, hence we can certify the generation of hybrid en¬ 


tangled entanglement. The corresponding experimental 
setup is shown in Fig. 2 (Tomo and Analysis sections). 
After a polarization analysis stage (two waveplates and 
a polarizer), each photon is sent to a g-plate and a sec¬ 
ond polarization analysis stage. In this configuration, the 
g-plate transfers the information initially written in the 
0AM subspace into the polarization state of the pho¬ 
ton that can be then analyzed with standard techniques 
[35, 36]. An hyper-complete set of measurements of po¬ 
larization and 0AM for both photons (overall 1296 set¬ 
tings) has been performed to fully reconstruct the den¬ 
sity matrix of the entangled photon pair in the sixteen¬ 
dimensional Hilbert space (2 dimensions for polarization 
and 2 for 0AM for each photon). Such characterization 
has been carried out for all the entangled pairs of differ¬ 
ent VV modes. The comparison between experimental 
and theoretical density matrices is reported in Fig. 3 A 
for two VV combinations. The complete set of experi¬ 
mental density matrices can be found in the Supplemen¬ 
tary Material. The quality of our states can be measured 
by calculating the fidelity F = 

between the experimental density matrix {pexp) and the 
corresponding theoretical one {ptheo)- The average fi¬ 
delity is: F = 0.97 ± 0.01. 

In order to quantify the intrasystem entanglement of 
VV states, we calculated the concurrence C for the single 
photon reduced density matrix in the polarization-0AM 
space after any projective measurement performed on the 
other photon. Figure 3B shows the concurrence distribu¬ 
tion for each photon when the other is projected over 34 
different states {xi}- As expected, the distribution is di¬ 
vided in two regions (blue and orange) corresponding to 
entangled ((7 = 1) and separable ((7 = 0) states respec¬ 
tively. Indeed when Xi contains a circular polarized state 
or a 0AM eigenstate the state of the other photon is sep¬ 
arable while in all the other cases it is a fully entangled 
state corresponding to a vector-vortex field. 

On the other hand, the intersystem entanglement be¬ 
tween two VV beams can be quantified by directly cal¬ 
culating the concurrence C for the corresponding den¬ 
sity matrix in the VV space. The average value for the 
concurrence over the five VV couples considered here is 
C = 0.91 ± 0.03 (see also Supplementary Materials). 

Beside the intriguing fundamental aspects of the en¬ 
tanglement, this feature is also a keystone in quantum 
based technologies. One of the most common ways to 
certify the entanglement is through a violation of a Bell 
inequality. Hence we performed a non-locality test for 
each pair of VV photons in order to show that entan¬ 
glement between complex vectorial fields constitutes an 
exploitable resource in quantum protocols. 

More in detail, we can map the bidimensional VV 
space of order m into a qubit encoding by defining 
|0) = l^m) and |1) = To test the entanglement 

we then measured the correlations between the two 
photons in the VV space. In a standard quantum 
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Figure 3. Experimental results. A) Experimental and theoretical density matrices for entangled VV beams of orders (mi = 
l,m 2 = 5) (on the left) and (mi = 3,m 2 = 5) (on the right). B) Polarization-0AM concurrence distributions. To quantify 
the intrasystem entanglement we calculated the concurrence distributions for the single photon reduced density matrix after 34 
different projective measurements performed on the other photon. All the distributions are divided in two regions corresponding 
to entangled (blue) and separable (orange) states. 


information scenario, two parties (Alice and Bob) 
perform measurements on one photon on the entangled 
pair respectively. According to this convention we 
performed a violation of the CHSH inequality [37] 
S = |-E(ao, bo) + E{ai^ bo) + E(ao^ bi) — -E(ai, Bi)\ < 2 
where cii{bj) are the outcomes of Alice’s (Bob’s) 
measurement settings i{j) and E{ai^bj) are the cor¬ 
relators between Alice and Bob measurements that 
quantify the probability of observing a coincidence 
in Alice and Bob detectors for a given measure¬ 
ment configuration. The maximum violation of the 
CHSH inequality allowed by quantum theory is given 
by S = 2\/2 corresponding to a maximally entan¬ 
gled state and to the following measurement bases: 
{|0),|1)} and ^{|0) + |1), |0) - |1)} for Alice and 
{cos(7r/8) |0) +sin( 7 r/ 8 ) |1), — sin( 7 r/ 8 ) |0) +cos( 7 r/ 8 ) |1)} 
and {sin(7r/8) |0) + cos( 7 r/ 8 ) |1), — cos( 7 r/ 8 ) |0) + 

sin( 7 r/ 8 ) |1)} for Bob. Being defined in the VV space, 
these measurement settings need the ability to perform 
projective measurements directly in the VV qubit 
Hilbert space. To this purpose, as mentioned above, 
we exploited a g-plate as an interface between VV 
and polarization spaces (Analysis stage in Fig. 2). 
For a given VV order m, a q-plate with q = m/2 
followed by a HWP and a polarizer (PBS) allows one 
to perform a projective measurement on a state of the 
form cos( 7 ) |0) + sin( 7 ) | 1 ), for any value of 7 . Different 


(mi, m 2 ) 

S (raw data) 

Scorr 

(1,1) 

(1.5) 
(1,10) 

(3.5) 
(1,-1) 

2.654 ± 0.009 

2.649 ± 0.013 

2.437 ± 0.010 

2.621 ± 0.016 

2.592 ± 0.010 

2.727 ± 0.009 

2.738 ± 0.014 

2.591 ± 0.011 

2.716 ± 0.017 

2.664 ± 0.010 


Table L Measured S parameter for CHSH inequalities using 
different orders of VV beams (first column) for raw data and 
for data corrected for dark counts (second and third column 
respectively). 

measurement settings for Alice and Bob correspond 
to different angles of their HWPs. The experimental 
values of S are reported in Table I and correspond to 
a violation of CHSH inequality by 39 to 73 standard 
deviations, without any correction due to dark counts. 

In this work we have fully investigated the hybrid en¬ 
tangled entanglement in a photonic system composed of 
two VV beams. Such a system shows indeed two differ¬ 
ent types of entanglement: an intrasystem entanglement 
between polarization and 0AM of each photon is respon¬ 
sible for the complex polarization pattern of VV beams, 
the two photons are also entangled with each other via 
an intersystem entanglement which lies at the basis of 
the non-locality concept. We investigated the structure 
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of our system by performing a full state tomography and 
quantified both types of entanglement by calculating the 
concurrence for the reduced density matrices in the sin¬ 
gle photon and VV modes Hilbert spaces respectively. 
Finally we have performed a non-locality test to prove 
that intersystem entanglement between entangled sub¬ 
system can be used as a resource in quantum protocols 
despite the complexity of the two subsystems involved. 
This study can pave the way towards a quantum enhance¬ 
ment in VV beams applications as well as the realization 
of quantum-information protocols that can take advan¬ 
tage of the combined action of both intrasystem and in¬ 
tersystem entanglement. 

We thank G. Rubino for her contribution in the prelim¬ 
inary stage of the experiment. This work was supported 
by PRIN (Programmi di ricerca di rilevante interesse 
nazionale) project AQUASIM and ERC-Starting Grant 
3D-QUEST (3D-Quantum Integrated Optical Simula¬ 
tion; Grant Agreement No. 307783): www.3dquest.eu. 
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EXPERIMENTAL APPARATUS INTERSYSTEM CONCURRENCE OF VV 

PHOTON PAIRS 


In the generation stage (Fig. 2 of main text) a 
pump beam of wavelength 397.5 nm is generated by a 
Tiisapphire mode-locked laser with a 76-MHz repetition 
rate that is upconverted by second harmonic generation 
(SHG). The pulses then impinge on a ^d-barium borate 
(BBO) crystal cut for type-II phase matching, generating 
pairs of entangled photons at 795 nm. Both photons are 
filtered in wavelength and spatial mode by using filters 
of A A = 3 nm and single-mode fibers, respectively. The 
unconverted pump beam is stopped just after the crys¬ 
tal. Figure 2 shows that a HWP can be inserted after 
one of the q-plates to flip the polarization handness in 
order to generate a 7r-mode: in particular, this was done 
for generating the VV pair (mi = 1, m 2 = — 1). The next 
step is the VV analysis that is performed for each pho¬ 
ton independently. Each photon is analyzed first in the 
polarization space by using a polarization-analysis stage 
(birefringent waveplates and a polarizing beam splitter) 
and then in the 0AM space. The 0AM analysis requires 
another q-plate working as an interface between polar¬ 
ization and 0AM. This q-plate is used to transfer the in¬ 
formation encoded in orbital angular momentum to the 
polarization degree of freedom that can be then easily 
analyzed using a second polarization analysis stage. Fi¬ 
nally, the photons are coupled into single mode fibers in 
order to detect only the states with null 0AM. 


We calculated the concurrence in the VV Hilbert space 
in order to quantify the intersystem entanglement be¬ 
tween the two photons. The results for the different pairs 
VV modes are shown in Table 1. 

These concurrence values were ob¬ 

tained after performing a projection over 
the four-dimensional subspace spanned by 

of the VV beam for the cases with m^ = 
1,3,5,10 and over the subspace spanned by 

> \rml} \^m2) > \Nl) \T^ti2) > \Nl) \T^m2)} 
the particular case of the 7r-mode (mi = l,m 2 = — 1). 


(mi, m 2 ) 

Concurrence 

(1,1) 

(1.5) 
(1,10) 

(3.5) 
(1,-1) 

0.949 ± 0.003 

0.906 ± 0.003 

0.863 ± 0.003 

0.908 ± 0.002 

0.914 ± 0.003 


TABLE 1. Intersystem concurrences for entangled VV beams 
of different orders 


QUANTUM STATE TOMOGRAPHY AND 
INTRASYSTEM CONCURRENCE 

Each of the entangled VV pairs has been analyzed 
via a full state tomography in the corresponding sixteen¬ 
dimensional polarization-0AM Hilbert space. In the fol¬ 
lowing pages are reported the experimental density ma¬ 
trices and the intrasystem concurrence distributions for 
all the entangled VV couples considered in this work. 
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FIG. 1. Experimental results. A) Experimental and theoretical density matrices for entangled vector vortex beams of order 
ml=l and m2=l. B) Concurrence distributions for the single photon reduced density matrix after 34 different projective 
measurements performed on the other photon. All the distributions are divided in two regions corresponding to entangled 
(blue) and separable (orange) states. 
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Entangled VV state: IV’1,5) 
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FIG. 2. Experimental results. A) Experimental and theoretical density matrices for entangled vector vortex beams of order 
ml=l and m2=5. B) Concurrence distributions for the single photon reduced density matrix after 34 different projective 
measurements performed on the other photon. All the distributions are divided in two regions corresponding to entangled 
(blue) and separable (orange) states. 
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Entangled VV state: IV’1,10) 
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FIG. 3. Experimental results. A) Experimental and theoretical density matrices for entangled vector vortex beams of order 
ml=l and m2=10. B) Concurrence distributions for the single photon reduced density matrix after 34 different projective 
measurements performed on the other photon. All the distributions are divided in two regions corresponding to entangled 
(blue) and separable (orange) states. 
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Entangled VV state: 1^3,5) = (l^s) li^s) + I’^s) l^s)) 

Re[p\ I'>T^[p\ 



Concurrence Concurrence 


FIG. 4. Experimental results. A) Experimental and theoretical density matrices for entangled vector vortex beams of order 
ml=3 and m2=5. B) Concurrence distributions for the single photon reduced density matrix after 34 different projective 
measurements performed on the other photon. All the distributions are divided in two regions corresponding to entangled 
(blue) and separable (orange) states. 
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FIG. 5. Experimental results. A) Experimental and theoretical density matrices for entangled vector vortex beams of order 
ml=l and m2=-l. B) Concurrence distributions for the single photon reduced density matrix after 34 different projective 
measurements performed on the other photon. All the distributions are divided in two regions corresponding to entangled 
(blue) and separable (orange) states. 

















